Satellite DNAs and the evolution of the multiple X1X2Y sex chromosomes in the wolf fish Hoplias malabaricus (Teleostei; Characiformes)

Multiple sex chromosomes usually arise from chromosomal rearrangements which involve ancestral sex chromosomes. There is a fundamental condition to be met for their long-term fixation: the meiosis must function, leading to the stability of the emerged system, mainly concerning the segregation of the sex multivalent. Here, we sought to analyze the degree of differentiation and meiotic pairing properties in the selected fish multiple sex chromosome system present in the wolf-fish Hoplias malabaricus (HMA). This species complex encompasses seven known karyotype forms (karyomorphs) where the karyomorph C (HMA-C) exhibits a nascent XY sex chromosomes from which the multiple X1X2Y system evolved in karyomorph HMA-D via a Y-autosome fusion. We combined genomic and cytogenetic approaches to analyze the satellite DNA (satDNA) content in the genome of HMA-D karyomorph and to investigate its potential contribution to X1X2Y sex chromosome differentiation. We revealed 56 satDNA monomers of which the majority was AT-rich and with repeat units longer than 100 bp. Seven out of 18 satDNA families chosen for chromosomal mapping by fluorescence in situ hybridization (FISH) formed detectable accumulation in at least one of the three sex chromosomes (X1, X2 and neo-Y). Nine satDNA monomers showed only two hybridization signals limited to HMA-D autosomes, and the two remaining ones provided no visible FISH signals. Out of seven satDNAs located on the HMA-D sex chromosomes, five mapped also to XY chromosomes of HMA-C. We showed that after the autosome-Y fusion event, the neo-Y chromosome has not substantially accumulated or eliminated satDNA sequences except for minor changes in the centromere-proximal region. Finally, based on the obtained FISHpatterns, we speculate on the possible contribution of satDNA to sex trivalent pairing and segregation.


Hybridization patterns of selected HmaSatDNA monomers in the HMA-D karyomorph
Among the 18 mapped HmaSatDNA monomers (marked by asterisk in Supplementary Table 1), only two (Hma-Sat14 and HmaSat56) did not produce specific signals on male and female mitotic chromosomes after FISH (data not shown).This outcome may be attributed to their potentially low abundance in the genome or their arrangement into clusters that are too small for effective cytological detection.The remaining examined satDNA monomers provided varied hybridization patterns (Figs. 1 and 2 and Supplementary Figs. 1 and 2).Only two hybridization signals, restricted to a certain autosome pair, were recorded for nine HmaSatDNA monomers, namely HmaSat11, HmaSat17, HmaSat25, HmaSat32, HmaSat34, HmaSat36, HmaSat42, HmaSat48, and Hma-Sat51, with HmaSat25 and HmaSat32 showing adjacent co-localized arrangement (Fig. 1 and Supplementary Figs. 1 and 2).On the other hand, FISH with HmaSat07 probe yielded strong (peri)centromeric signals exclusively on sex chromosomes, extended towards the p-and q-arms of the X 1 and neo-Y, respectively (Fig. 2 and Supplementary Figs. 1 and 2).HmaSat07 was hence the only sex chromosome-specific satDNA monomer revealed in the present study.Finally, the remaining six satDNA monomers displayed FISH signals on multiple autosomes, neo-Y, and either X 1 or X 2 sex chromosomes.HmaSat01 and HmaSat18 exhibited highly similar patterns, with strong clusters being detected in the telomeric regions of many autosomes along with X 2 and neo-Y chromosomes.On the other hand, HmaSat03, HmaSat04 and HmaSat05 shared the pattern of (peri)centromeric signals, being present on many (HmaSat04) or almost all (HmaSat03 and HmaSat05) chromosomes of the complement, including neo-Y and either X 2 (HmaSat03) or X 1 (HmaSat04 and HmaSat05) chromosome.A specific pattern was produced by the probe for HmaSat02, with prominent clusters at many terminal and interstitial regions across the entire chromosome set (Fig. 2).

Patterns of HmaSatDNA monomers on the meiotic sex trivalent formed by X 1 X 2 Y multiple sex chromosomes of HMA-D
To better understand the organization of HmaSatDNA monomers on the X 1 X 2 Y multiple sex chromosomes of HMA-D during their meiotic pairing into the trivalent structure, we sought to map the seven HmaSatDNAs that display detectable clusters on these sex chromosomes (Fig. 2).A prominent cluster was found for the Hma-Sat01 and HmaSat18 monomers, corresponding to the segment formed by the terminal regions of the X 2 and neo-Y chromosomes (Fig. 3).Four sites that resided on the X 1 , X 2 , and neo-Y chromosomes were displayed by the HmaSat02-specific probe.In contrast, the HmaSat03 probe displayed two signals, one of which populated the q-arms of the neo-Y chromosome and the second located in the centromeric region of the X 2 chromosome (Fig. 3).All the three remaining HmaSatDNA monomers (HmaSat04, HmaSat05, and HmaSat07) presented two FISH signals in the (peri)centromeric region of the X 1 and neo-Y chromosomes (Fig. 3).

Chromosomal location of HmaSatDNA monomers in HMA-C
We mapped the subset of seven HmaSatDNA monomers, which hybridized to HMA-D sex chromosomes, also to metaphases of H. malabaricus males of karyomorph C (HMA-C) to examine their possible presence on nascent XY sex chromosomes (Fig. 4 and Supplementary Fig. 3).Most of FISH signals resembled the patterns revealed on the HMA-D sex chromosomes, with the exception of HmaSat01 (Fig. 2).The same satDNA monomers found on the X 1 and the (peri)centromeric region of the neo-Y, (HmaSat02; HmaSat04; HmaSat05; HmaSat07) were located in the (peri)centromeric regions of the corresponding X and Y chromosomes in the karyomorph HMA-C (Fig. 4 and Supplementary Fig. 3).In a similar way, the HmaSat01 and HmaSat02 monomers, with signals on the X 2 and neo-Y (HmaSat01), or all multiple sex chromosome (HmaSat02) of the HMA-D, also formed telomeric signals in both nascent sex chromosomes of the HMA-C (Fig. 4 and Supplementary Fig. 3).In turn,  the HmaSat03 and HmaSat18 monomers were the only ones without detectable hybridization signals on the X and Y chromosomes of HMA-C.In the case of HmaSat03, the situation particularly reflects the fact that clusters of this satDNA reside on the X 2 sex chromosome of HMA-D, i.e. the newly added chromosomal part (Fig. 2), therefore this satDNA joined the sex chromosome system only after the Y-autosome fusion.On the other hand, HmaSat18 was entirely absent from the HMA-C chromosomes, at least under the FISH resolution (Fig. 4 and Supplementary Fig. 3).Lastly, it is worth noting, that the probe corresponding to HmaSat05 monomer marked the centromeric regions of all but one chromosome pair in the HMA-C complement, which, after comparison with the HMA-D pattern (particularly the lack of signal on X 2 ), indicates that the sole chromosome pair without the signal represents highly probably the linkage group (pair 20) which has been incorporated into the X 1 X 2 Y system in karyomorph HMA-D.

Minimum spanning trees (MST)
We selected two satDNA monomers to generate minimum spanning trees: (1) the HmaSat02, as the satDNA with the most pronounced hybridization pattern, with telomeric and interstitial signals located on all chromosomes of the HMA-D complement, including multiple sex chromosome system, and (2) the HmaSat18, which displayed terminal signals on several autosomes and in the X 2 and neo-Y chromosomes of this karyomorph.The MST of HmaSat02 revealed the predominance of three abundant haplotypes, besides the presence of many monomers with few mutations.All of them displayed a very similar abundance to all haplotypes identified in males and females (Fig. 5).The MST of HmaSat18 revealed the presence of many haplotypes.Although the predominant ones were shared between males and females, unique male-specific haplotypes were observed (Fig. 6).

Discussion
In this work, we studied the satDNA content and its chromosomal distribution in one cytotype (karyomorph D) of the wolf-fish H. malabaricus with multiple sex chromosome system of the X 1 X 2 Y type.The main purposes of this study were to investigate, (1) whether certain satDNAs could be associated with the breakpoint regions involved in the Y-autosome fusion which gave rise to the mentioned multiple sex chromosome system, and (2) whether the neo-Y chromosome already underwent a significant degree of differentiation that could be reflected by the satDNA accumulation.The opportunity to compare the satDNA distribution also with the patterns on the XY sex chromosomes of closely related karyomorph (HMA-C), from which X 1 X 2 Y of HMA-D originated via Y-autosome fusion 74 , qualifies H. malabaricus as a useful study model for these research questions.
A recent burst in satellitome studies has led to significant advancements in our understanding of the teleost satDNA evolutionary dynamics and its possible impact on various genomic processes 66,[68][69][70]77 . Thu far, a particular attention has been given to characiform fishes 64,[66][67][68]77 .Among the characiform satellitomes, the diversity of satDNA ranges from 30 (Piaractus mesopotamicus) to an exceptional diversity of 164 monomers (Megaleporinus macrocephalus) 68,78 .In the present study, we revealed 56 satDNA monomers in the genome of H. malabaricus-karyomorph D (Supplementary Table 1), which falls within the above-mentioned range reported for characiforms.We further revealed an extraordinary wide range of satDNA monomer sizes, ranging from 6 to 2774 bp, with the predominance of those exhibiting higher proportion of AT-pairs (> 50%) (Supplementary Table 1), similarly to what has been found in some other vertebrates such as mammals and amphibians [79][80][81] , but differing from reptiles and birds with a predominantly CG-rich content 82,83 . Sixteen outof 18 chromosomally mapped HmaSatDNA monomers showed detectable clusters pointing on their organization into high-copy number arrays.We observed the following three major distribution patterns: (1) only two hybridization signals restricted to a certain autosome pair (nine monomers; Fig. 1); (2) clusters located on the X 1 , X 2 , and/or neo-Y sex chromosomes in addition to several-to-all autosomes (six monomers; Figs. 2, 7), and (3) a single satDNA, HmaSat07, accumulated exclusively in the pericentromeric regions and surrounding area on the X 1 and neo-Y sex chromosome (Figs.2, 7).
Intriguingly, although the most abundant satDNA sequence in the genome very often occupies centromeric regions of all chromosomes, implying its functional significance for this crucial chromosomal region 84 , our data on distribution patterns imply this is not the case for HmaSat01, the most abundant satDNA monomer in HMA-D.Previously, the 5SHindIII-DNA satDNA, corresponding to the herein characterized HmaSat03 monomer, has been proposed to play a role in centromere function or structure in this fish 76 , however, our current analysis showed that other satDNA, namely HmaSat05, is a more likely candidate for this role based on the higher proportion of (peri)centromeric signals per chromosome set (Fig. 2).Yet, this monomer does not form detectable clusters on each centromeric region in the complement, with a single autosome pair (putatively the autosomal addition creating X 1 X 2 Y system in HMA-D) lacking the signal in HMA-C, and more chromosomes without the signal in HMA-D (Figs. 1, 2, 4).Our present findings therefore suggest a gradual change in the (peri) centromeric satDNA content in H. malabaricus, with analogous pattern being recently reported for African Nothobranchius killifishes 85 .In H. malabaricus, this process seems to be considerably fast given that HmaSat05 monomer occupies the centromeric regions only in HMA-C and HMA-D, while it is not present in these regions on chromosomes of closely related karyomorphs A, F, and G (Supplementary Fig. 4).It is generally thought that the mechanisms driving fast turnover of centromeric repeats and associated proteins may be responsible for a reproductive barrier between diverging populations of the same or closely related species, due to hybrid incompatibility related to chromosome segregation [86][87][88] .This, together with possible incompatibilities related to different sex chromosome systems 89 , might have contributed to the previously reported absence of hybrids between H. malabaricus karyomorphs 90 , thus supporting the proposal that the karyomorphs of H. malabaricus represent already distinct evolutionary units 71 .A functional role of HmaSat05 monomer in HMA-C and HMA-D centromeres will need to be further assessed by testing the interaction between this satDNA and key centromeric proteins such as CENP-A 52,58 .
From this point on, we will attempt to address the two main queries raised beforehand, starting with the inquiry of whether changes in satDNAs content accompanied the rearrangement(s) involved in the origin of the X 1 X 2 Y multiple sex chromosome system present in HMA-D.
The association between repetitive DNA accumulation and the occurrence of chromosomal rearrangements along with overall elevated genome instability has been extensively recorded in diverse lineages across the tree of life 91,92 , including humans 93 .In the case of X 1 X 2 Y sex chromosome system of HMA-D, most neo-Y-linked clusters of HmaSatDNA monomers are concentrated in the centromeric region and its surrounding area (Fig. 7).From the herein mapped monomers, only HmaSat02 and HmaSat03 repeats populate the putative breakpoint junction of   the Y-autosome tandem fusion (Figs. 2 and 3), which has been inferred from previous WCP experiments 72 , and comparisons of the position of repetitive DNA clusters between the X 2 and neo-Y (HMA-D) and the ancestral XY sex chromosomes in HMA-C 74 .Whether or not these sequences directly contributed to the rearrangement remains an open question.The other monomers, HmaSat04 and HmaSat05, forming strong (peri)centromeric accumulations on X 1 and neo-Y, are likely absent from the ancestor of the X 2 chromosome as demonstrated by the lack of FISH signals in Figs. 2 and 3. Contrary to other studies 61,62 , which evidenced a high rate of elimination and turnover of satDNAs populating neo-sex chromosomes, our findings indicate that both HmaSat04 and HmaSat05 satDNAs were not excised during the Y-autosome fusion event or lost through the subsequent process of centromere inactivation/elimination (reviewed by Barra and Fachinetti 94 , and McKinley and Cheeseman 95) .
The second question was whether or not the neo-Y chromosome underwent a significant degree of differentiation, which could be inferred from the level of satDNA accumulation.For this, we took the opportunity of direct comparison of satDNA landscapes between the multiple X 11 X 2 Y sex chromosomes of HMA-D and standard XY sex chromosome system of HMA-C from which the system of HMA-D evolved.Species with heteromorphic sex chromosomes usually exhibit sex-linked variations in the accumulation of certain repetitive DNA sequences; they can be either substantially enriched or even limited to sex chromosome(s) 31,32,43,61,62 , which has been reported also in teleosts 64,67,68 .When narrowed down to multiple sex chromosomes, those so far revealed in teleost fishes show, in the overwhelming majority, low degree of repeat accumulation 2 .In the case of X 1 X 2 Y system of HMA-D, our results corroborate the previous findings from other FISH-based studies [73][74][75] by showing non-elevated degree of satDNA accumulation, with a limited difference between the synteny blocks of X 1 /X 2 versus neo-Y chromosome (Supplementary Table 1; Fig. 7).Minor changes linked to neo-Y satDNA content may be related to the loss of the HmaSat01 accumulation from the terminal portion of q-arms (on both neo-Y and X 1 chromosome), based on the comparison with patterns found in the XY sex chromosome system of HMA-C (see Fig. 7).Nevertheless, it is not possible to rule out the alternative hypothesis that HmaSat01 repeats populated/amplified to cytologically detectable extent on the XY sex chromosomes of HMA-C more recently, i.e. the sex chromosome-autosome fusion in HMA-D might predate this accumulation.
The causal relationship between repetitive DNA accumulation and recombination suppression on sex chromosomes has been a matter of a long-standing debate.It seems more likely that repetitive DNA accumulates in the region of suppressed recombination and then it contributes to further expansion of this region, to genetic decay, and to sex chromosome heteromorphy 4,[30][31][32][33][34][35]96 . Imprtant question is then: what happens during the first meiotic division when sex chromosomes need to pair efficiently despite possible differences in size, morphology and/or repetitive DNA/heterochromatin content?Accumulation of repeated DNAs and heterochromatin could hamper the establishment and stability of the trivalent, preventing thereby its proper segregation, leading to high frequency of unbalanced gametes and thus being avoided by natural selection [20][21][22]28,97,98 . Bearing his in mind, we speculate that the lack of major differences in satDNA accumulation between the members of X 1 X 2 Y sex chromosome system of HMA-D, along with generally low degree of repeat accumulation might have contributed to faithful pairing and segregation of these sex chromosomes (Fig. 7).The presence of sex chromosome-specific repeat segments in the area of assumed sequence differentiation, both in the XY sex chromosome system of HMA-C and X 1 X 2 Y system of HMA-D, might facilitate smooth sex chromosome pairing within this region, without asynapsis, as has been already proposed for other studied teleost fishes 99,100 .

Conclusions
Taken together, our present study provides novel insight into the aspects linked to emergence and differentiation of representative fish multiple sex chromosome system-X 1 X 2 Y present in H. malabaricus, karyomorph D. We showed that the satDNA content remained rather unchanged during the evolution of this system after Y-autosome fusion.Last but not least, by unravelling the satellitome features of this H. malabaricus karyomorph, we provided new tools for further studies on satDNA evolution in characiform fishes.

Sampling and chromosome preparation
We analyzed 20 individuals of H. malabaricus karyomorph HMA-C (11 males; 9 females) and 17 individuals of karyomorph HMA-D (10 males; 7 females).HMA-C individuals were collected in Bento Gomes River (Poconé, Mato Grosso State, Brazil) and HMA-D in the Monjolinho stream (São Carlos, São Paulo State, Brazil).The sex of sampled individuals was identified by morphological analysis of gonadal tissues, and each individual was assigned to specific karyomorph based on the characteristic cytogenetic profiles (e.g., 2n, karyotype structure) uncovered in previous studies 71,73,74 .Sampling was authorized by the Brazilian Environmental Agency ICMBIO/ SISBIO (License 48628-14) and SISGEN (A96FF09).Mitotic chromosomes were obtained from kidney biopsies following the protocol by Bertollo et al. 101 , while meiotic spreads were prepared from testes of HMA-D males according to Kligerman and Bloom 102 .Briefly, we treated the animals with benzocaine and applied a colchicine solution (0.025%) for 40 min.Subsequently, we excised kidney samples and testicular tissue, which were fragmented and treated with a hypotonic solution (0.075 M KCl) for 20 min at 37 °C.Both mitotic cell suspensions and fragmented testes were preserved in Carnoy I fixative (3:1 methanol/glacial acetic acid).Chromosome slides were stained with a 10% Giemsa solution (pH 6.8).All experiments followed the ethical guidelines sanctioned by the Ethics Committee on Animal Experimentation of the Universidade Federal de São Carlos (Process Number 7994170423).The authors complied with ARRIVE guidelines.

DNA extraction and next generation sequencing
Genomic DNA (gDNA) of one male and one female of HMA-D was extracted using liver tissue by the standard phenol-chloroform-isoamylalcohol procedure 103  www.nature.com/scientificreports/prepared from the isolated HMA-D gDNAs and sequenced on the BGISEQ-500 platform at BGI (BGI Shenzhen Corporation, Shenzhen, China).The sequencing yielded a total of 1.55 Gb and 1.4 Gb of raw reads for males and females, respectively.The genomic reads obtained were deposited in the Sequence Read Archive (SRA) under accession numbers SRR27786732 (male) and SRR27786731 (female).

Analysis of H. malabaricus satDNAs
We characterized the satellite DNA content of HMA-D male and female independently.First, we examined the read quality (Q > 20 for all nucleotides), by using the Trimmomatic software version 3.0 104 .Then, a subset of 2 × 500,000 reads has been randomly subsampled and usedas input for the TAREAN analysis 51 .All recovered satDNAs were removed from the original library by using the DeconSeq software version 0.4.3 105 , and from the remaining dataset, a new randomly selected subset of 2 × 500,000 reads was created to pass through a second round of the TAREAN analysis.These steps (iterations) were repeated until no novel satDNAs were found.

Estimating the abundance and genetic distances of satDNAs
The abundance of each satDNA (proportion of reads containing the same repeat vs. the entire pool of reads) was estimated by selecting 2 × 500,000 reads from each genomic library (HMA-D female and HMA-D male).Then, using the "cross-match" tool in the RepeatMasker version 4.1.6106 and a custom python script (https:// github.com/ fjrui zruano/ ngspr otoco ls/ blob/ master/ repeat_ masker_ run_ big.py, accessed on 20/01/2023), the reads were aligned against each respective catalog.We used the script calcDivergenceFromAlign.py in the RepeatMasker 106 to calculate the genetic distance of each satDNA (Kimura-2-parameter).All satDNAs were named "HmaSat" and numbered according to their decreasing abundance in the genome.Finally, to check for the possible presence of male-biased satDNAs, we calculated the male/female (M/F) abundance ratio of each HmaSat monomer given by RepeatMasker and selected those with an M/F ratio greater than 1.0.

Minimum spanning trees (MST)
To achieve the accurate scores of monomer diversity, we used PHYLOVIZ version 2.0 107 to construct the minimum spanning trees (MSTs) for two satDNA monomers which were found on the sex chromosomes of HMA-D and, at the same time, their monomer sizes fell into the required range for the analysis, i.e. below 150 pb (reflecting the size of the sequencing raw reads): HmaSat02 and HmaSat18.First, we subsampled 2 × 5,000,000 reads from each genomic library (males and females HMA-D) and aligned the reads using Bowtie2 version 1.3.1 108 against the two aforementioned satDNA monomers.Then, we removed the central area that corresponded to one monomer and, in order to avoid sequencing error bias by CD-HIT version 4.8.1 109 , we also removed singletons (monomers detected only once) and aligned each satDNA independently using MUSCLE version 5.0 110 .

Primer design and DNA amplification via polymerase chain reaction (PCR)
Out of the 56 HmaSatDNA monomers found by our analysis (see Supplementary Table 1), we selected 18 representative monomers (five most abundant ones and those with some differences in the abundance between sexes; indicated by asterisk in Supplementary Table 1) for physical mapping.For this purpose, we manually designed primers for each consensus monomer, using both multiple primer analyzer (Thermofisher) and Oligo-Calc (Biotools) to check for the hairpin formation and self-annealing.The PCR conditions followed Kretschmer et al. 67 and it was performed in a T100 Bio-Rad Thermal Cycler (Bio-Rad Laboratories).Briefly, the amplification reaction consisted of an initial denaturation at 95 °C for 7 min, followed by 34 cycles with denaturation at 95 °C (45 s), annealing temperatures varying from 50 to 62 °C (60 s), extension at 72 °C (60 s), and the final extension at 72 °C (7 min).All PCR products were quantified using a ThermoFisher NanoDrop spectrophotometer (ThermoFisher Scientific) and submitted to electrophoresis with a 1% or 2% agarose gel (Supplementary Fig. 5).

Fluorescence in situ hybridization (FISH)
To examine the chromosomal location of HmaSatDNA monomers, with the emphasis on their landscapes on the nascent XX/XY (HMA-C) and multiple X 1 X 1 X 2 X 2 /X 1 X 2 Y sex chromosomes (HMA-D), we performed FISH on both male and female metaphase spreads of the two karyomorphs.To unambiguously identify the sex chromosomes in both karyomorphs, we took the advantage of the previously presence of 18S ribosomal DNA (rDNA) on the neo-Y and X 1 chromosome of HMA-D, and XY chromosomes of HMA-C 74 and the characteristic distribution patterns of one satDNA monomer 76 .In addition, the X 2 chromosome was identified according to its morphology and size, as it is the smallest submetacentric chromosome in the karyotype 71,74 .For FISH, we used the 18S rDNA probe (1400 pb-long segment of the 18S rRNA gene from H. malabaricus) was previously isolated by Cioffi et al. 111 .Finally, we also used the sex chromosome-specific marker, HmaSat07 satDNA described in the present study, in complementary experiments.Results from the sequential experiments (i.e. on the same metaphases) with 18S rDNA or HmaSat07, and each individual satDNA are provided in Supplementary Figs. 6 and 7.The HmaSat DNAs that displayed signals on the mitotic X 1 , X 2 , or neo-Y chromosomes of HMA-D, were mapped also onto male meiotic spreads in prophase I of this karyomorph and further onto HMA-C metaphases.All the 18 tested satDNA monomers along with the 18S rDNA fragment were labeled with Atto-488-dUTP; Atto-425-dUTP (green fluorescence), or Atto-550-dUTP (red fluorescence), using a nick-translation kit (Jena Bioscience, Jena, Germany) according to the manufacturer's instructions.The hybridization procedures were conducted

Fig. 3 .
Fig. 3. Meiotic chromosomes of male Hoplias malabaricus karyomorph D (HMA-D) after FISH with a set of HmaSatDNA monomers located both on autosomes and sex chromosomes.Arrows indicate the sex trivalent.Bar = 5 μm.

Fig. 4 .
Fig. 4. Male mitotic metaphases of Hoplias malabaricus, karyomorph C (HMA-C) after FISH with a set of HmaSatDNA monomers that produced detectable signals on X 1 X 2 Y sex chromosomes of karyomorph D (HMA-D).SatDNA family names are indicated on the right bottom, in green (Atto-488-dUTP) or red (Atto-550-dUTP).The arrows point on the X and Y chromosomes.Bar = 5 μm.

Fig. 5 .
Fig. 5. Linear MST of HmaSat02 satDNA obtained from female (red) and male (green) reads.The diameter of the circle is proportional to the abundance of the haplotype.Black circles represent one base of divergence between haplotypes.

Fig. 6 .
Fig. 6.Linear MST of HmaSat18 satDNA obtained from female (red) and male (green) reads.The diameter of the circle is proportional to the abundance of the haplotype.Black circles represent one base of divergence between haplotypes.